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Abstract. Evidence for strong electron-phonon coupling in high-T c cuprates is 
reviewed, with emphasis on the electron and phonon spectral functions. Effects 
due to the interplay between the Coulomb and electron-phonon interactions are 
studied. For weakly doped cuprates, the phonon self-energy is strongly reduced due to 
correlation effects, while there is no corresponding strong reduction for the electron self- 
energy. Polaron formation is studied, focusing on effects of Coulomb interaction and 
antiferromagnetic correlations. It is argued that experimental indications of polaron 
formation in undoped cuprates are due to a strong electron-phonon interaction for 
these systems. 
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1. Introduction 

High-T c cuprates show a large number of interesting features, apart from the 
exceptionally large superconducting transition temperature T c . They exhibit 
antiferromagnetic pQ, pseudogap [2J, marginal Fermi liquid [3] and ordinary Fermi 
liquid phases in addition to the superconducting phase. After the high-T c cuprates 
had been discovered [3], there was initially much interest in the electron-phonon 
interaction (EPI). It was, however, soon concluded that the EPI is too weak to 
explain superconductivity alone, in particular <i-wave superconductivity, and the interest 
focused on purely electronic models of these compounds. More recently, there has been 
substantial experimental evidence that the EPI plays an appreciable role for a number of 
properties. Certain phonons show a large softening and broadening under doping [HE], 
suggesting a strong interaction with doped holes. This is, for instance, seen for the 
so-called half-breathing copper-oxygen bond stretching phonon, apical oxygen phonons 
and the oxygen Bi g buckling phonon. Photoemission spectroscopy (PES) experiments 
show the formation of small polarons for the undoped cuprates [7j, and a kink in the 
nodal k-direction also suggests strong EPI [8]. While there is only a weak isotope 
effect on T c for optimally doped samples, a strong isotope effect has been seen away 
from optimum doping [9]. Recent STM work suggests that a phonon mode plays 
a role in superconductivity [TO], although other interpretations are possible |TT]. In 
particular, an isotope effect has been observed [10]. While the phonon contribution to 
superconductivity remains unclear, it seems clear that phonons can be important for 
other properties. 

The EPI has been studied very extensively in the local density approximation 
(LDA) [12] of the density functional formalism [13], which is particularly appropriate for 
systems where correlation effects are not very strong. This approach has been shown 
to be very successful for conventional superconductors [131 IH2, ETUI H7j . For cuprates 
[T8] a rather weak EPI was found, which alone would not be sufficient to explain the 
superconductivity [19]. However, the calculated width [20] of the half-breathing phonon 
is an order of magnitude smaller than the reported experimental value [21] , raising some 
questions about the accuracy of the LDA in this context [22J . This is one of the reasons 
that the interest has recently focused on whether the interplay between the Coulomb 
interaction and the EPI can explain experimental signs of a strong EPI. 

Due to the important effects of the Coulomb interaction in these systems, models 
such as the Hubbard and t-J models are often used. In these models important phonons 
couple to charge fluctuations. Since charge fluctuations are strongly suppressed in the 
cuprates by the Coulomb interaction, an important issue is if this could mean that the 
EPI is actually suppressed. We discuss this issue extensively below. 

In the so-called sudden approximation, angular resolved photoemission spec- 
troscopy (ARPES) can be directly related to the one-electron Green's function. If 
superconductivity is due to bosons coupling to electrons and forming electron pairs, 
this coupling should show up in the one-electron Green's function. Due to the high 
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energy- and k-resolution that can now be obtained in ARPES, a lot of interest has 
focused on ARPES recently, and we address these issues below. 

ARPES experiments strongly indicate that small polarons are formed for undoped 
cuprates and there are signs of strong phonon side bands [TJ. This indicates that 
there is a strong EPI for these systems. For weakly underdoped or optimally doped 
cuprates, ARPES experiments show quasiparticles, suggesting that there are no small 
polarons formed in these cases. However, there is still substantial spectral weight in 
the energy range where phonon side bands would be expected, suggesting that the 
EPI is still substantial. There has been extensive work on polarons and bipolarons 
in metals, treating both electronic properties in general and superconductivity [23J. 
Since experiment suggest that small polarons are not formed at dopings relevant for 
superconductivity, we here focus on polaron formation for insulating systems. 

Due to the great interest in cuprates, there have been many reviews covering 
different aspects of these systems [H El El El El El HH [2H [251 EH [271 EH [29l [301 EU 

H21E21E2!. 

2. Models 

2.1. Coulomb interaction and hopping 

The Coulomb interaction plays an important role in the cuprates. A frequently used 
model for describing this is the three-band model [35], which includes a Cu x 2 — y 2 3d 
orbital and two O orbitals in a Cu02 plane. The model includes the Cu-0 hopping 
integrals and the Coulomb interaction between two electrons on the Cu orbital. 



where and So are the energies of the Cu and O atoms, respectively. 5 describes the 
O atom positions in the unit cell and runs over (a/2, 0) and (0, a/2) in the second term 
and over (±a/2,0) and (0, ±a/2) in the third term, where a is the lattice parameter. 
P_a = -P s , P s = 1 for 5 = (a/2,0) and P s = -1 for 5 = (0,a/2). c\ a creates a Cu 
electron in cell i with spin a, a\ Sa creates an O electron and n a = c* ia c ia . U is the 
Coulomb interaction and t p d is a hopping integral. 

From this model the t-J model can be derived [36j, where each site corresponds to 
a Cu atom in the Cu0 2 plane. In the undoped system, corresponding to all Cu atoms 
being in d 9 configurations, each site is occupied by one hole. In a hole doped system, 
the holes go primarily onto the O sites. Such an O hole forms a Zhang- Rice singlet with 
a Cu hole [36] • A Zhang- Rice singlet is described by an empty site in the t-J model. 
The corresponding Hamiltonian is 




(1) 



+ t pd j2 p s( 



4a a i6v + H - C + U 
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where (ij) refers to a sum over nearest neighbor pairs, and c\ a creates a spin a hole on 
site i if this site was previously empty. Sj is the spin and = J2 a ^L^o- i s the number 
of holes on site i. 

Alternatively, the one-band Hubbard model is often used 

#Hub = -t ^{c\ a c ja + H.c.) + U ^ntT^u- (3) 

(ij)cr i 

The t-J model can also be derived from the Hubbard model in the large U limit if 
certain terms are neglected [37] . 

2.2. Electron-phonon interaction 
We introduce the Hamiltonian for a set of phonons 

q 

where 6 q annihilates a phonon with the frequency cj q and the wave vector q 
we write the coupling to the phonons as 

#el-ph = -L 9 ^ ^i+q/k^q + & -q)> 

kq 

where N is the number of cells. 

2.2.1. Holstein phonons Often, the electron-phonon interaction is treated in a Holstein 
model, where there is an on-site coupling to one local Einstein phonon per site. This 
corresponds to a k and q independent coupling 

0Hoi(k, q) = 9o, (6) 
where go is the coupling constant, and a q independent phonon frequency cc^Hoi = oj^. 

2.2.2. Breathing phonons The breathing (oxygen bond-stretching) phonons have 
attracted much interest due to the observation of an anomalous softening and broadening 
of these phonons when the system is doped [3] . That these phonons may have a strong 
coupling can be understood by noticing that the formation of the Zhang- Rice singlet in 
the t-J model involves a large energy of the order of several eV. For a system without 
phonons and a fixed number of doped holes, this energy only enters as an uninteresting 
constant. If phonons are added, however, the singlet energy can be modulated by the 
phonons. This is the case for the breathing phonons, where the O atoms in the Cu02 
plane move in the direction of the Cu atoms, thereby changing the bond lengths (see 
Fig. CD). This directly modulates the Cu-0 hopping integrals t p d (in a three-band model) 
determining the Zhang-Rice singlet energy and leads to a substantial coupling. This has 
been discussed by several groups [38j ESI El HH H2J 152 HI] . A general formula for this 
coupling was given in Ref. j31] . considering both the modulation of the Cu-0 hopping 
integrals and shifts of the levels due to Coulomb interactions. It was found that the main 



(4) 

. Generally 

(5) 



Interplay between electron-phonon and Coulomb interactions in cuprates 



5 




Figure 1. Half-breathing (upper left), breathing (lower left), Ai g (upper middle), Bi g 
(lower middle) and apical oxygen Of (right) phonon modes. In the left and middle 
figures, the small filled circles show Cu atoms and the large circles O atoms in the 
CuC>2 plane. In the right figure [5j, the large black spheres show apical O atoms. 

coupling is an on-site coupling due to the modulation of the Cu-0 hopping integrals. 
One reason for this result is that the hopping integrals in the t-J model, obtained after 
the O levels have been projected out, are about an order of magnitude smaller than 
the on-site singlet energy. This strongly favors the on-site electron-phonon interaction 
over the coupling to the t-J hopping integrals [15]. Below we therefore neglect the off- 
site interaction. If we furthermore assume that the vibration of the Cu atom can be 
neglected due to its larger mass, the oxygen phonon eigenvectors can be approximated 
as e° = sin(ag Q ,/2)-^/sin 2 (ag :r ./2) + sin 2 (ag y /2), where a = x or y, describing the motions 
in the directions of the nearest neighbor Cu atoms. Then the coupling becomes 

g-Br(k, q) ~ yj sin 2 (ag«/2) + sih^ag^), (7) 
i.e., the coupling depends on q but not on k in this approximation [46J. 

2.2.3. Ai g and B\ g phonons There has been substantial interest in a B lff phonon at 
about 42 meV [4TI l4~8l |4~9] , involving a motion perpendicular to the Cu02 plane of the 
O atoms in this plane. These atoms can move out of phase, leading to a B l9 phonon, or 
in phase, leading to an Ai g phonon (see Fig. [1]). Devereaux and coworkers [50l [5TJ [52] 
studied these phonons in the three-band model. They considered the case when there 
is an electrical field, E z , perpendicular to the Cu02 plane, due to different valencies of 
the ions on the two sides of a plane. Since the O atoms move perpendicular to the plane 
for the Ai 5 and B lg phonons, these phonons couple to such a field. For the case of only 
nearest neighbor hopping, this leads to the coupling [52J 

9A lg ,B lg (k, q) ~ E z [B x( f>l(k + q)0 x (k) + BM + q)<k,(k)], (8) 
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where </>. = T^(k)/V# 2 (k) + 2 (k), E(k) = \e d - e |/2 + v^d - ^o) 2 /4 + ^ 2 (k), 
f2 2 (k) = £ 2 (k) + ty(k) and t a (k) = 2tsin(a/c a /2). Using a spring model, Devereaux et 
al. [53J obtained the eigenvectors [B x ,B y ] = [cos(aq y /2), — cos(ag x ./2)]/M(q) for the 
Bi 9 mode, where M(q) = ^ cos 2 (aq y / 2) + cos 2 (aq x /2). Putting the Cu mass equal to 
infinity, their model gives [B x , B y ] = [cos(aq x /2), cos(aq y /2)]/M(q) for the Ai fl mode. 

In contrast to the breathing phonons, these modes have a strong k dependence. 
For q = 0, the coupling due to the B lg phonon is entirely off-site, while the A lg phonon 
has a substantial on-site coupling. 

A different approach was taken by Jepsen et al. [31]. They studied a six-band 
model of the LDA band structure and focused on the coupling due to the modulation of 
hopping integrals. This coupling is zero for a completely flat CuC>2 plane but becomes 
finite for a dimpled plane. They obtained a coupling 

, , a(k x + q x ) ak x a(k y + q y ) ak y . . 

9B lg (k, q) ~ cos cos — ~ cos ^ — ^cos^. (9) 

This coupling tends to emphasize small values of |k| and |k + q| more then Eq. (jSJ), due 
to the appearance of cos-functions rather than sin-functions in t a (k) entering in Eq. (jH]). 



2.2-4- Apical phonons Neutron scattering experiments show that several apical oxygen 
phonons (see Fig. CP broaden and soften when a cuprate is doped [3J. This coupling 
has been calculated for La 2 Cu0 4 using a shell model [55] . Due to the ionicity of the 
O atoms, the electrostatic part of the coupling is expected to be particularly strong. 
This is in particular true for the undoped system, which is an insulator, leading to a 
poor screening of the electrostatic interaction. It is therefore important not to perform 
this calculation using the LDA, since LDA gives a metallic system and too efficient 
screening. Due to the electrostatic nature of the coupling, and due to small q vectors 
playing an important role, the coupling is expected to be dominated by the on-site part 
of the coupling, i.e., 

9 Apical (k, 

q)~0A(q)- (io) 

Calculations showed that this coupling is indeed rather strong |55j . 



3. Weak coupling and noninteracting electrons 

3.1. Electron self-energy 

The electron-phonon interaction is often studied assuming that the electrons are 
noninteracting. This is a quite unrealistic assumption for the cuprates, where the 
electron-electron interaction is crucial. Below, we nevertheless describe some of the 
results [561 E3 EH] f° r this case, since they provide a basis for discussing similarities 
and deviations for strongly correlated systems. The electrons are described by the 
Hamiltonian 

H non = ^2e k cl a c ka , (11) 

kcr 
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Figure 2. The spectral function (left), calculated for different values of the bare energy 
£k (in units of w p h), and the dispersion (right), comparing EDC and MDC results with 
the bare dispersion. The inset in the left hgure shows the real and imaginary part of 
£. All results are obtained for the model in Eqs. (|6lllj) and they have been given a 
Lorentzian broadening (FWHM) of 2d. The parameters are A = 0.5, u> p h/B = 0.05 
and = 0.05. 



where is the energy for the wave vector k and a is a spin index. The electrons are 
assumed to couple to the phonons via the Holstein model (Q. We calculate the retarded 
electron self-energy to lowest order in the coupling g. For T = 0, it is given by [56J 

S ( k > <") = ^9 2 Et =7 + W(£k+q) J , (12) 

N u + uj ph - e k+q + id u-uj ph - e k+q + id 

where N is the number of sites, f(e) is the Fermi function and 6 is a positive infinitesimal 
(later small) quantity. We assume that N(e) = 1/B is constant, where N(e) is the 
density of states (DOS) per spin and B is the band width. The band is assumed to be 
half-filled and to extend from -B/2 to B/2. Then the sum in Eq. (|T2l) can be performed 
exactly, giving 

w , \ x x ri uj + uj ph + B/2 + i5 u-u ph + i5 

s k > w) = -Xu ph In y - — + In p , 13 

2 uj + u; ph + id u — u ph — B/2 + id 

where 

A = -^JV(O). (14) 
E is k-independent in this approximation. Assuming that o; p h <C 5/2, we obtain 

I I Pin 

0, if w ph < |w| < B/2. 

Results for S are shown in the inset of Fig. [2j At T = 0, a hole with a binding 
energy larger than the phonon frequency cj p h can excite a phonon and simultaneously 
be scattered. |ImS| is therefore large for u < — co> p h- If the system is (s-wave) 
superconducting, scattering is only possible for binding energies larger than tUph + A, 
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where A is the gap, since the scattered electron has at least the binding energy A 
Related to the abrupt onset of Im E there is a logarithmic singularity in Re S. The left 
part of Fig. [2] shows the spectral function 

A(v n 1 |ImS(k,o;)| 

1 ' ' 7r[^-£ k -ReS(k,^)] 2 + [ImS(k,^)] 2 ' 1 ' 

For simplicity, in Fig. [2] we have assumed a linear dispersion over the energy range of 
interest, e k = vp(k — kp), where t> F and kp are the Fermi velocity and wave vector, 
respectively. A(k,uj) shows a narrow peak at E k = £ k /(l + ^) if l-^kl ^ ^ph and it has 
a broad peak at E k = e k if |e k | 3> uj ph . The electron-phonon coupling therefore leads to 
a change of slope 

\dE k /dk\ Ek =o 1 



(17) 



|d£; k /(ik]| Bk | >a;ph 1 + A' 
We also define the quasiparticle strength 

z(k) = i-(9s(k,^)/^)U £k = TTa' (18) 

where the second equality is only valid for |i? k | <C uj ph . 



3.2. Energy and momentum distribution curves 

Figure [2] shows energy distribution curves (EDC), where A(k, uj) is shown as a function 
of uj for a fixed value of k. For each value of k we can determine the uj for which A(\c, uj) 
has its maximum. From this we obtain a dispersion relation E^ shown by the red curve 
(EDC) in the right part of the figure. From the results for A(k, uj) it is clear that the 
EDC dispersion jumps from E^/uj^ ~ —1 for e^/uj^ « —2 to E-^/uj^ w —2.1 for 
£k = /^ P h ~ —2.3. The EDC dispersion is reduced by a factor (1 + A) close to E^ = 0, 
as discussed below Eq. (fl6j) . and the dispersion is further reduced as _E k approaches 

Alternatively, we can study momentum distribution curves (MDC), showing A(k, uj) 
as a function of k for a fixed uj. In particular, we can find the k for which A(k,u) 
has its maximum (uj fixed). This MDC dispersion relation is shown in the right part 
of Fig. [2j As an illustration, the vertical lines in the left part of the figure show 
the energies uj/uj^ = —1 and -1.25. The maximum (among the k values shown in 
Fig. [2]) of A(k, uj/ujph = —1) is obtained for e^/ujp^ = —2, while the maximum of 
A(k, uj/uj ph = —1.25) is obtained for e^/uj-p^ = —1.55. Over a certain frequency range, 
an increase in \lj\ then leads to a decrease in |e k |. The result is the S-like shape of the 
MDC dispersion around uj/lj p ^ = — 1. Well away from this energy, the EDC and MDC 
curves agree for this k-independent self-energy. 

A(k, uj) tends to show a dip at uj ~ — co> p h, in particular if > £ k > — uj p ^ [58]. This 
is the combined effect of a logarithmic singularity in Re E(k, uj) at uj — — co» p h, which 
makes the first term in the denominator in Eq. ( fl6|) large, and the sudden onset of Im 
E(k, uj) below uj = — Wph, which makes the numerator large. As discussed above, the 
singularity is shifted to uj p ^ + A in a (s-wave) superconductor. At the same time the 
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effect becomes stronger, since weight is piled up at the onset of Im E(k, u), leading 
to a more pronounced structure in Re S(k, u). As discussed by Sandvik et al. |59j, 
this effect is particularly pronounced for an s-wave superconductor and less strong for 
a (i-wave superconductor. 

Similarly, the phonon self-energy Il(q, uj) can be calculated to lowest order in 
g, where q is the phonon wave vector. Im Il(q, uj) is of particular interest, since 
it determines the phonon width. For a system of noninteracting electrons and a 
nondegenerate phonon mode, Allen [60], El] showed that 

2 Im n(c; ph ) = 2nh 2 cul h N(0)\, (19) 

where 2 Im Il(q, Wph) is the full width at half maximum (FWHM) and iV(0) is the 
density of states (DOS) per spin at the Fermi energy. Knowledge of the FWHM of the 
phonons then gives a possibility of estimating A. 



3.3. Electron-phonon coupling in the 2d Holstein model 

Two-dimensional (2d) correlated models with EPI are often compared with the 2d 
Holstein model to determine the effects of correlation on the EPI. A 2d Holstein model 
at half-filling with only nearest neighbor hopping is unstable to an infinitesimal EPI 
due to perfect nesting. Therefore the comparison is often made to a Holstein model 
with just a single electron at the bottom of the band [62j [63] . Often a t-J model doped 
with one hole is studied, suggesting similarities with a Holstein model with a single 
electron. The half-filled Holstein model, however, is of particular interest, since the 
relevant antibonding Cu-0 band in the cuprates is close to half-filling. 



3.3.1. Weak- coupling limit We consider the limit g <C t and uj p h <C t for a nearest 
neighbor hopping t with a single electron Holstein model. Using the self-energy in 
Eq. ( TT2T) , we can define an effective electron-phonon coupling A = Ao via the quasiparticle 
weight Z [Eq. (TT8]) ] or the effective mass 



d 2 E k /dk 2 



d 2 S\n/ rfk 2 



= TTT' (20) 

k=0 1 + A o 



where 



Ao = -r^—. (21) 

Both methods give the same A for this model in the weak-coupling limit [64J . Instead of 
the Holstein model with a single electron, we can study the half-filled model assuming a 
constant density of states (DOS) and calculating S(k, uj) according to Eq. ( TT31) . Defining 
A via the expression for Z(k) recovers the A defined in Eq. ( THI) and leads to A = n\ , 
since iV(0) = l/(8t). The increase in A is partly due to the fact that the self-energy at 
the bottom of the band only has contributions from higher states while in the half-filled 
case there are contributions from both higher and lower states and partly due to the 
DOS at the bottom of the band being smaller than the average DOS |65j . 
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3.3.2. Polaron formation Polaron formation has been studied extensively for the 
Holstein model [661 E3 E3 [691 EH EH [631 E2]- For noninteracting electrons, polaron 
formation is often associated with bipolaron formation [7H [72] . Since on-site bipolaron 
formation is strongly suppressed by the large on-site U relevant for cuprates, we focus 
on polaron formation here. We use Z — * as the criterion for polaron formation. For 
the nearest neighbor Holstein model with a single electron and = 0.0125iy this was 
found to happen for A c = 1.2 [63J, where W is the band width and A is here defined 
as A = 2g 2 /W, corresponding to the assumption N(0) = 1/W in Eq. ( 1T4"1) . For the 
half-filled case with a semi-elliptical DOS and using the dynamical mean-field theory 
(DMFT) [73] is was found that A c = 0.33 [65J, again using cUph = 0.0125W 7 and putting 
N(0) = 1/W in the definition of A. Similar results were found by several other groups 
[69~1 1701 17T1 172] . As in the weak coupling limit, there is a large difference between the 
single electron and half-filled cases [65] . 

To understand this difference, we consider polaron formation in the adiabatic limit 
by comparing free electron states with states of perfectly localized electrons [7U 165] . 
The energy for free electrons is Ef ree = —At per electron for the single electron case and 
Ef r ee = —16t/(37r) ~ —1.72 for the half-filled case. In the localized case, E\ oc = —g 2 /ujp\ 1 
per electron for both cases. For simplicity, we assume that polarons form when 
l-Eiod > | -Efree I- This leads to a large A c = 1 for the single electron case and a much 
smaller A c = 0.42 for the half-filled case, in rather good agreement with more accurate 
calculations. The large difference between the two cases is due to the large difference in 
kinetic energy per electron. In the single electron case, the electron is at the bottom of 
the band and has the maximum (absolute) kinetic energy, while in the half-filled case 
the average kinetic energy is much smaller. The electron-phonon interaction energy can 
then win more easily and lead to polaron formation. 

4. Experimental results 

4-1. Phonon softening and width 

There have been extensive studies of phonons in high-T c cuprates using neutron 
scattering. For reviews, see, e.g., Pintschovius and Reichardt [5], Pintschovius [6], 
Egami and Billinge [25] and Egami [31]. Here we focus on results of particular interest 
for the EPI and in particular La2_ I Sr :r Cu04, where the most complete results have been 
obtained. Figure [3] shows results of Chaplot et al. [75] for undoped (left hand side) and 
doped (right hand side) L^-xSr^CuC^. The solid curves have been obtained within 
a shell model [75]. The shell model calculations give a rather accurate description of 
almost all phonon branches. 

One striking exception is the highest mode of Ai symmetry for the doped system. 
This is the so-called half-breathing phonon, which is a bond-stretching vibration of the 
oxygen atoms in the Cu02 plane, as shown in Fig.[TJ This phonon is rather well described 
by the shell model for the undoped system. However, doping leads to a strong softening 
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Figure 3. Phonon dispersion curves for La2Cu04 (left hand side) and Lai.gSro.iCu04 
(right hand side). The solid curves were calculated in a shell model by Chaplot et 
al. [75] . The dash-dotted lines were obtained after inclusion of a quadrupolar force 
constant. The figure illustrates that in the doped compound (right hand side) the 
half-breathing mode in the Ai symmetry at 20 THz and several apical oxygen modes 
in the Ai symmetry at 12-16 THz are softened. 

half way along the (1,0,0) direction. This softening is anomalous in the sense that it 
is not captured by the shell model. Anomalies of such bond-stretching phonons have 
been observed by several groups for both high-T c cuprates [76| [TTj [78| [2^ [79j [8(7t [8Tt [82] 
and other compounds [83j [MJ [85]. The anomalous behaviour of this phonon is also 
illustrated by the large broadening in the doped system, as is shown in Fig. [U Thus 
Pintschovius and Braden [21] found the FWHM intrinsic broadening for q = (0.5, 0, 0) to 
be 1.2 THz=5 meV for 15% doping. It is unlikely that this width is caused by disorder, 
since the width of the q = (0.5, 0.5, 0) quadrupolar mode (where two of O atoms around 
a Cu move towards Cu and two O atoms move away from Cu) is resolution limited [86J. 
Both the softening and the width indicate that this phonon couples strongly to doped 
holes. Using the formula [Eq. (fl9|) ] of Allen [601 16T] for noninteracting electrons and the 
density of states iV(0) ~ 0.66 states per eV and spin [87], the electron-phonon coupling 
can be estimated to be A ~ 0.1 — 0.3 for the (half-)breathing phonons. In view of the 
arguments in Sec. 15.11 this estimate and similar estimates below should be taken with 
caution. 

Reznik et al. [88] found a strong anomaly for the bond-stretching phonon around 
q = (0.25, 0, 0) for several cuprates. This anomaly was particularly large for systems 
with static stripe order but it was also seen for systems where stripe order has not 
been observed, e.g., superconducting La 2 _ 2; Sr :E Cu04. This was interpreted in terms of 
a coupling to a dynamic charge inhomogeneity. 

The Of phonon of A x symmetry at about 17 THz in the undoped system and with 
the reduced wave vector (001) shows strong softening under doping (see Fig. [3]). As 
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Figure 4. Width of the oxygen bond-stretching phonon for La2- x Sr x Cu04 as function 
of the reduced wave vector (£,0,0) and for doping x = (left), x = 0.1 (middle) 
and x = 0.15,0.30 (right). The full line shows the experimental resolution including 
focusing effects. The figure illustrates the large broadening for the doped system and 
the small intrinsic broadening for the undoped (x = 0) system (after Pintschovius [6]). 



shown to the right in Fig. [IJ this phonon mainly involves the movement of apical O. 
The softening of this phonon was predicted by Falter et al. [89], [90] before being observed 
experimentally. The width of the Of phonon is 16 meV [5] E], suggesting a very strong 
coupling to this mode. 

There has also been a substantial interest in a Bi 9 phonon involving out-of-plane and 
out-of-phase vibrations of oxygen atoms in the C11O2 plane. This phonon has an energy 
of about 42 meV for YE^CusC^^. The phonon has been studied using both Raman 
[HI [921 ESI S3 EH EI] and neutron scattering [33, EH]- This phonon shows an interesting 
change of frequency and width as the compound is cooled below T c [91] [92] [93] HH] |4"9] . 
and a Fano line shape is observed in Raman scattering [92] [93]. Fitting of the line 
shape and changes of the spectrum as the sample is cooled, leads to estimates of the 
electron-phonon coupling in the range A ~ 0.02 — 0.06 [37J MM EI] f° r YBa 2 Cu 3 07_5, 
comparable to theoretical estimates A = 0.02 from band structure calculations |94j in 
the local density approximation. Qualitatively similar but smaller effects have also been 
observed for other phonons [3H1 EH] . 

4-2. Photoemission spectra 

Photoemission spectra have been extensively reviewed [26] [29] [32] [33] , and here we only 
show a few typical results of particular interest in the context of the EPI. 

Much interest in the EPI has been created by the observation of a kink in the 
experimental electron dispersion for several cuprates [96] [97] [98] [99] 11001 [8] 11011 1102j . 
Some typical results are shown in Fig. E] Lanzara et al. [8] emphasized that such a 
kink is found for three different families of compounds (LSCO, Bi2201 and Bi2212), for 
different dopings and both below and above T c . Other groups obtained similar results 
for Bi compounds, but disagreed about whether there is a kink above T c |101[ 1102] or not 
[991 1100] . The structure for LSCO, both below and above T c , is more pronounced than 
for B2212 (see Fig. E]). For noninteracting electrons, the ratio of the slopes below and 
above the kink is expected to be given by the dimensionless electron-phonon coupling 




Figure 5. Dispersion in the (0, 0)-(7r, tt) direction for La2-5SraCu04 (LSCO), 
Bi 2 Sr 2 Cu0 6 (Bi2201) and Bi 2 Sr 2 CaCu 2 8 (Bi2212) for different dopings S and 
temperatures T. The results in (a) and (b) were obtained for T = 20 K and in 
(c) for T = 30 K. The red arrow shows the energy of the q — (tt, 0) half-breathing 
phonon and the black arrows the energies of the kinks. Panel (f) shows the change 
1 + A' of the slope at the kink (after Lanzara et al. [8]). 



1 + A (Eq. (|17p ). Figure [5f shows the change of slope 1 + A'. Here A' may be different from 
A, since even the states at large binding energy may not show the "bare" dispersion. If 
the theory for noninteracting electrons were applicable (Eq. (TlTI) ). this would suggest a 
coupling of the order of A ~ 1. 

While the early measurements showed only one kink at about 70 meV, later work 
found several structures at smaller binding energies [1031 1104] - Estimates of Re S(k, u) 
were extracted from experiment by assuming that the underlying "bare" dispersion is of 
second order in |k| — k-p over the range of interest. The second derivative of Re £(k, to) 
for La2_ a: Sr a: Cu04 then has structures at about 40-46 and 58-63 meV and possibly at 
23-29 and 75-85 meV, suggesting that there is coupling to bosons at these energies. 

The dispersion of the quasiparticle for Bi 2 Sr 2 Ca n _iCu n 02n+4 away from the nodal 
direction and, in particular, in the antinodal (tt, 0) direction has been extensively studied 
[1051 M, UM UM HDSl \M\. For Bi2212 and Bi2223 below T c there is a very strong 
structure away from the nodal direction. As discussed in Sec. EJ if the kink is due to a 
coupling to a mode at energy ub, the kink is expected to appear at roughly u B + A in 
the superconducting state [58j I108j . where A is the gap. This suggests ujb ~ 40 meV 
for optimally doped Bi2212. The boson could be the Bi ff buckling phonon mode [107j . 
which has roughly the right energy. Alternatively, the coupling [106] could be to the 
so-called resonance peak, seen in inelastic neutron scattering in the superconducting 
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Energy (eV) k (O.OHtc.jc) 

Figure 6. a) The k = (tt/2, 7r/2) spectrum of undoped Ca2Cu02Cl2 as well as fits to 
a Gaussian (solid red curve) and a Lorentzian (dashed blue curve). The maximum of 
the broad feature is denoted by A and its onset by B. b) Dispersion of A and B along 
the nodal direction as well as the different values of the chemical potential fj, for a large 
number of samples (after Shen et al. [7]). 



states of Bi2212 and YBCO pm EOl HH EE]. Recently, a broad structure has also 
been seen at 40-70 meV the inelastic neutron scattering spectrum for L^-zSr^CuCU 

4-3. Chemical potential. Polarons 

A long-standing problem has been the position of the chemical potential in undoped 
or strongly underdoped cuprates as well as the interpretation of the PES spectra for 
these systems. Figure shows the spectrum of undoped Ca2Cu02Cl2 at the top of the 
band [k = (7r/2, 7r/2)]. The spectrum has a broad feature centered at A and an onset 
at B. Figure [6b shows the dispersion of these features. The dispersion of A (Fig. [6b) 
agrees well with the extended t-J model [26]. It has therefore often been assumed 
that A represents a quasiparticle which is very strongly broadened by some unknown 
mechanism. 

Shen et al. [7J pointed out, however, that the peak shape is not Lorentzian, as 
would be expected from a life-time broadening, but Gaussian. Even more seriously, 
a quasiparticle at the top of the valence band of an insulator cannot decay into an 
electron- hole pair, and one would expect the width to be small [7J. In Fig. [6b, each 
horizontal line shows the chemical potential of a specific sample. Since the system is an 
insulator, one might expect the chemical potential to be pinned to impurities or defects 
and therefore depend on sample preparation. The figure shows, however, that although 
the value of the chemical potential is highly sample dependent, it is always at least about 
0.45 eV above the peak A. Shen et al. [7] pointed out that all these puzzling features 
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can be explained by assuming polaronic behaviour due to strong coupling to bosons. 
Peak A is then a boson side-band, explaining its large width and Gaussian shape, and B 
represents the quasiparticle, having too small a weight to be seen experimentally. This 
explains why the chemical potential is never lower than B in Fig. [6b, since this is the top 
of the valence band in the new interpretation. It has furthermore been observed that 
the width of the broad peak increases substantially with T, providing further support 
to the interpretation in terms of coupling to bosons [114] . 

As the doping is increased, a quasiparticle is observed in Ca 2 - :c Na a ;Cu02Cl2 for 
x = 0.1 [7J. The quasiparticle weight is small, however, and there is still much weight in 
the energy range of the side band. This suggests that the effects of the electron-boson 
interaction remain important for the doped system, although (small) polaronic effects 
are not seen anymore. For La2- x Sr x Cu04 this is found to happen at smaller dopings 
(x = 0.03) [115] . but it has been proposed that this could be due to phase separation 

4-4- Isotope effects 

Isotope effects provide a good indication for electron-phonon interactions. In particular, 
the interest has focused on the superconductivity transition temperature T c . A review of 
this work is given by Franck [DJ. Generally, a small oxygen isotope effect was found for 
compounds having the composition which gives the highest T c in that family. However, 
the isotope effect is often much larger for systems where T c is suppressed by some 
substitutions. 

Khasanov et al. |117j observed an isotope effect in the penetration depth of nearly 
optimally doped YBa2Cu 3 7 _5, using a muon-spin rotation technique, which allows a 
direct observation of the penetration. These results were interpreted in terms of an 
isotope effect in the carrier mass due to a strong electron-phonon interaction violating 
Migdal's theorem |117j . 

Gweon et al. |118j studied the oxygen isotope effect on the PES spectra for 
optimally doped Bi 2 Sr2CaCu20 8+ 5 and found a very large effect, involving shifts as 
large as 30-40 meV. Douglas et al. |119j in contrast found no large isotope effect, but 
concluded that their measurements are not inconsistent with a conventional isotope shift 
of the order of 3 meV. 

4-5. Scanning tunneling spectroscopy 

There has recently been much work based on scanning tunneling spectroscopy (STM), 
observing strong spatial modulations and a checkerboard structure [1201 11211 11221 [TO] . Of 
particular interest here is the work of Lee et al. [ID] , studying Bi 2 Sr 2 CaCu208+,5. From 
the second derivative of the tunneling current d 2 I/dV 2 they determined the spatially 
dependent gap A(r) and the energy of a bosonic mode Q(r) appearing in the spectrum. 
The mode energy showed the isotope effect expected for phonons involving primarily 
O atoms. Based on this and on the doping independence of O(r), it was concluded 
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that the mode is an O phonon. Its average energy was found to be 52 meV, with a 
substantial spatial variation. They found anticorrelation between fl(r) and A(r). They 
concluded that the results could be due to i) a heterogeneity of the frequency and 
coupling constants of a pairing- related phonon causing a disorder in A(r), ii) inelastic 
scattering of tunneling electrons by phonons unrelated to superconductivity [H] or iii) 
a competing phase coupling to phonons and causing the anticorrelation between f2(r) 
and A(r). 



5. Interplay between Coulomb and electron-phonon interactions 

In systems where the Coulomb repulsion U is important, the effects of the electron- 
phonon interaction can be strongly influenced. This can easily be seen for the Holstein- 
Hubbard model [Eqs. (J3J, ED] - The Hamiltonian can be transformed so that the phonons 
couple to the deviation of the site occupancies from their average. If U is small, 
the number of electrons on a given site fluctuates strongly, and there are substantial 
deviations from the average, even for a half-filled system. However, as U is increased, 
the fluctuations are reduced. For U large and close to half-filling, most sites then 
have exactly one electron. For these sites there is no electron-phonon coupling. This 
suggests that the effects of the electron-phonon coupling are reduced as U is increased. 
The problem is, however, substantially more complicated, as discussed below. For the 
cuprates, phonon frequencies are typically smaller than electronic energies. We therefore 
do not discuss the antiadiabtic limit, where phonon frequencies are much larger than 
electronic energies and quite different effects can be found [1531 [72] . 



5.1. Sum rules 

We first discuss the interplay between Coulomb and electron-phonon interactions in 
terms of sum-rules for the imaginary parts of the electron and phonon self-energies. 
We consider the t-J model [Eq. ([2])] and the electron-phonon coupling in Eq. §5§ for k 
independent coupling #(k, q) = g(q), which results in an on-site coupling. 



5.1.1. Phonon self-energy We first consider the phonon self-energy Il(q, to). This can 
be written as 

n(q ' u) - iTwm^m^y (22) 

where x(q, u) is the charge-charge response function and Do(c\,uj) is the free phonon 
Green's function. Khaliullin and Horsch [39j showed that this function satisfies a sum 
rule (at T = 0) 

i r°° 

— ]T/ \lm X (q,uj)\duj = 25(l~5)N, (23) 

where 5 is the doping. This result is suppressed by a factor 25(1 — 5) compared with the 
result for noninteracting electrons in a half-filled band. Since x(q, a;) becomes small for 
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small 5, the denominator in Eq. (|2"2"j) is not very important, and the sum-rule in Eq. (T2"3l 
also applies approximately to Il(q, uj)/g^ [123J 

i i r°° 

— / |ImII(q,u;)|du;«25(l-5), (24) 

7riV q^O 9 1 J -°° 

Since typically 5 ~ 0.1, this implies that the softening and width of a phonon due to the 
creation of electron-hole pairs is drastically reduced. A formula of Allen [60, 61J is often 
used to estimate the electron-phonon coupling A from the phonon width 2 Im II(q, u). 
This formula is derived for noninteracting electrons and it neglects the strong reduction 
of n in Eq. ( 1241) . Its use for high-T c cuprates may therefore substantially underestimate 
A. 



5.1.2. Electron self-energy To derive a similar sum rule for the electron self-energy, we 
define a Green's function 

G (M)= ,_ £k ! k s(M) . (25) 

where is a weight and E(k, z) is the electron self-energy. The z-independent part of 
S is included in the energy e^, so that S(k, z) ~ b^/ z for large z. By studying the large 
z behaviour of G(k, z) and S(k, z), Rosch and Gunnarsson [123J related a sum rule for 
E(k, z) to moments of the electron-phonon interaction part of the Hamiltonian. For the 
t-J model and 5 = 0, the sum rule takes a very simple form 

1 f° 1 

- / ImE ep (k,^ - i0 + )du = ^J2 l^i 2 = 9*> ( 26 ) 

^ J — oo 1 * q 

where S ep is the difference in self-energy for the system with and without EPI. Since 
the inverse photoemission (IPES) spectrum has zero weight for 5 = 0, the integration 
extends only to u = 0. This result is shown in Fig. [Th. For noninteracting electrons and 
to lowest order (in g^), a similar sum rule for |ImS ep (k, uj — i0 + )\ (averaging over k) gives 
the right hand side g 2 /2 for integration up to uj = or g 2 for integration over all energies 
(see Fig. [7h). Thus the effect of the electron-phonon coupling is strongly reduced by the 
small doping for the phonon self-energy but not for the electron self-energy, at least not 
in a sum rule sense. 

The sum-rule in Eq. ( 126]) can be understood as follows. The electron Green's 
function describes the creation of a hole (or an electron). The phonons couple strongly 
to the charge of this hole, even if charge fluctuations are strongly suppressed elsewhere in 
the system. On the other hand, if a phonon is created, there is only a small fraction 5 of 
singlets that can respond. As a result, the electron-phonon interaction can be expected 
to appear to be a factor of l/(c<5) stronger in Re S(q, a>) than in Im H(q,u) (phonon 
width), where c ~ 2 — 4 depends on the assumptions about the cu-dependencies of Im E 
and Im IT. 

Although the arguments above show that the sum-rule (J26l) should not go to zero 
for 5 — > 0, the result is, nevertheless, nontrivial. The derivation depends on the coupling 
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Figure 7. a) Weights obtained by integrating Im £ op (k, u> — i0 + ) / tt over the indicated 
frequency intervals for the half-filled Hubbard and undoped t-J models. Also shown 
are the results for the Hubbard model without vertex corrections (no VC) and the 
lowest-order (k-averaged) result for the U = Hubbard model (non-int.). For the 
t-J model, the photoemission spectrum has been shifted by -U/2 and for the U = 
Hubbard model the PES and inverse PES (IPES) spectra have been shifted by -U/2 
and U/2, respectively, b) U dependence of the absolute value of the total and partial 
sum rules for Im E op using DMFT with A = 0.0025. The dotted lines indicate the 
expected small-C/ behaviour. W is the U = band width (after Rosch et al. |124| ) . 
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Figure 8. Contribution to S ep . The full and dashed lines are electrons and phonon 
Green's functions, respectively, and the circle is a vertex correction. 



being on-site and we have found no simple result for an off-site coupling. The right hand 
side is independent of k, t and J and it remains proportional to g 2 even for large g. 

Similar results can be derived for the half-filled Hubbard model in the large U limit 
[124] . The PES spectrum is expected to be close to the result for the t-J model, but the 
sum rule differs by a factor of two, due to the different integrated weights of the total 
(PES and IPES) spectra of the two models (see Fig.[7h). Integrating Im S ep (k, u—i0 + )/n 
over all frequencies for the large U Hubbard model gives the sum rule g 2 . Therefore 
there is a contribution — 3g 2 close to uo = 0, showing how the EPI reduces a large positive 
contribution at uj = already present for g = 0. Going from U — to large U, the 
sum rules for the PES or IPES part increase by a factor of four. This is illustrated 
in Fig. [7b- Interestingly, a substantial part of this change has already happened when 
U is comparable to the U = band width W [124] . The phonon sum rule is further 
illustrated in Sec. 16.11 




2468 10 2468 10 

U U 

Figure 9. Real part of g(p, q) as a function of U for (a) p = (— tt, 0) and (b) 
p = (— tt/2, tt/2). The value of q is indicated by the shape of the symbol. The solid 
circles are Monte Carlo results and the open symbols show results from perturbation 
theory (after Huang et al. [128] ). 



5.2. Vertex corrections 

In a diagrammatic language, an important contribution to S ep is shown in Fig. [HI 
expressed in terms of a vertex correction V. Figure [7^ shows that the neglect of V 
(no VC) in the half-filled large U Hubbard model leads to a large violation of the sum 
rules in the previous section. It is therefore interesting to study vertex corrections 
[I251II2EllI271H2glH2giH301H5I] . 

In particular, Kulic and Zeyher [1251 I127j have argued that in the large U limit, 
vertex corrections favor small angle scattering over back scattering. This would have the 
important effect of favoring c?-wave superconductivity. Huang et al. |128j performed a 
determinantal quantum Monte-Carlo (QMC) calculation [132] for the Holstein-Hubbard 
model for a 8 x 8 cluster for T = 0.5, where the hopping integral was put t — 1 and the 
band width is W = 8. The occupancy was n = 0.88. They defined an effective coupling 

g(p,q) = T(p,q) v /Z(p)Z(p + q), (27) 

where T(p, q) is the vertex correction, p and p + q stand for the momenta and imaginary 
frequencies of the incoming and outgoing electrons, respectively, and Z is the quasi- 
particle weight. The frequencies were put at their minimum value 7rT. Results are 
shown in Fig. It illustrates how g(p,q) is reduced relative to its U = value for 
all parameters shown, due to the reduction of the quasiparticle weights Z. Small angle 
scattering [q = (tt/4, tt/4)], however, is favored over back scattering [q = (jr, tt)] for large 
values of U, which helps <i-wave superconductivity. Some questions have been raised 
related to the large value of T which has to be used in the QMC calculations [1291 D-30J. 

5.3. Effects due to antiferromagnetic correlations 

Equation ( J261) shows that Im S ep (k, uj) is not suppressed by correlation effects in a 
sum rule sense. The sum rule does not tell us, however, how the contributions are 
distributed in frequency. Correlation effects tend to reduce the dispersion and provide 
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more low-lying excitations to which a quasi-particle could couple. This tends to increase 
Re |£ ep (k, u)\. On the other hand, spectral weight is removed from the quasi-particles 
and shifted away from the Fermi energy by correlation effects. This tends to have the 
opposite effect. These two effects need not cancel, as discussed below. 

5.3.1. Exact diagonalization There have been a large number of studies of the 
Holstein-Hubbard, Holstein-t-J or related models based on exact diagonalization 
[1551 HM [1351 [1361 H571 11381 [1391 [HDl M, ED 032] . Here we show some typical results. 
Zhong and Schiittler [134] studied a Holstein-Hubbard model i/Hub + -^Hoi (Eqs. EJ [6]) in 
the adiabatic limit (u; p h = 0) using exact diagonalization for an 8-site (a/8 x a/8) cluster 
with one doped hole, i.e., 7 electrons. They found that the system goes from a delocalized 
state to a polaronic state for A ~ 0.2 — 0.4 in the range U/t ~ 8 — 12. This was compared 
with a calculation for a spin-polarized system with one hole, where the antiferromagnetic 
spin correlations are removed. In this case the transition to a polaronic state occurred 
for A ~ 0.96. This was interpreted in terms of antiferromagnetic correlations strongly 
reducing the coupling A needed to obtain polaronic behaviour. Similar conclusions were 
obtained for larger cluster and using small but finite phonon frequencies by Fehske et 
al. [137[ I138[ D-40J, who studied Holstein-t-J and Holstein-Hubbard clusters. These 
calculations showed that the adiabatic approximation overestimates the tendency to 
(small) polaron formation. The quasiparticle weight and dispersion were shown to 
be reduced as the coupling was increased. Bauml et al. [140] studied the optical 
conductivity and argued that the midinfrared peak in cuprates is mainly of electronic 
origin while it is mainly of polaronic origin in nickelates. 

5.3.2. Self- consistent Born approximation Ramsak et al. [62J used a self-consistent 
Born approximation (SCBA) [1431 11441 11451 11461 1147] for a Holstein-t- J model, where 
a hole is assumed to couple to magnons and phonons. The problem is treated in a 
diagrammatic approach, neglecting crossing phonon or magnon lines. The electron- 
phonon part of the self-energy was written as 

S cl _ ph (k, u) = — ^2 glG(k -q,u- u ph ), (28) 
q 

where N is the number of sites, g q a coupling constant, G(k,u) the electron Green's 
function and u; p h the phonon frequency. Equation ([281 assumes that the electron-phonon 
coupling is not very strong and it neglects an indirect contribution, AE m [148] . via the 
coupling of magnons to the changes of the Green's function due to the electron-phonon 
coupling. To interpret their calculation, Ramsak et al. [62] argued that the main 
contribution to £ c i- P h is due to the coherent part G co h(k, ui) = Z(k)/(u — e^) of the 
Green's function, where Z(k) is the quasi-particle weight. They furthermore assumed 
that the quasi-particle energy in the absence of phonons can be parameterized in 
terms of effective masses m\\ and m_|_, i.e., ~ = k?J2m\\ + k"j_/2m±, where fey and 
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k± are measured from (tt/2,tt/2) in the (tt/2,tt/2) — > (0,0) and (7r/2,7r/2) — > (vr, 0) 
directions, respectively. This gives 

1 /• _ Z(k-q)o 2 
S cl _ ph (k, w = 4-— / d V 1 4 7 , (29) 

where the factor of four is due to the presence of four hole pockets at (±7r/2, ±7r/2) and 
the g-integrations over the surroundings of a hole pocket have been extended to infinity. 
This leads to a simple approximate formula for the enhancement of the mass m* in the 
(tt/2, 7r/2) — > (0,0) direction due to the electron-phonon interaction, 



— = (1 - 4A Z 2 ^^)- 1 = (1 + A), (30) 
my mo 

where Z is approximated by Z(ir/2,ir/2), A [Eq. ( 12T1) ] is the A obtained for a 
Holstein model with just one electron at the bottom of the band and mo = l/|2t| 
is the corresponding mass. Equation ( 1301) contains an extra factor Z |148j . obtained 
from solving the Dyson equation for the quasi-particle energy with electron-phonon 
interaction. As shown in Fig.dOl Eq. ( 1301) then agrees rather well with the full calculation 
for a large range of J/t values. The figure illustrates that the use of the coherent part 
[Eq. fl29l) ] is a rather good approximation for large J/t, while for small J/t it substantially 
underestimates A, primarily due to the neglect of AE m . 

The expression for m* /mil in Eq. (1301 is reduced by Z 2 but enhanced by the large 
effective mass factor y/m\\m±/mo, representing the reduced mobility of the carriers due 
to antiferromagnetic correlations. Thus the antiferromagnetic correlations greatly help 
polaron formation [62]. For large J/t, the net effect is an increase of A/Ao- For J/t = 0.3, 
the factor Z 2 y / m]jmX/mo is actually smaller than unity, but A/Ao is still enhanced, due 
to the factor of four in Eq. (130]) . resulting from the four hole pockets. 

5.3.3. Quantum Monte-Carlo calculations Mishchenko and Nagaosa [63] studied 
an undoped infinite Holstein-t-J model using a diagrammatic Monte-Carlo method. 
In contrast to the treatment above, they included diagrams with crossing phonon 
propagators and only neglected diagrams with magnon propagators crossed by phonon 
or other magnon propagators. Using the parameters u^/t = 0.1 and J/t = 0.3, they 
found self-trapping for A ~ 1.2 for an electron at the bottom of the band in the Holstein 
model but already for A ~ 0.4 in the Holstein-t-J model [63]. It was concluded that 
the antiferromagnetic ground-state and the coupling to magnons help the formation of 
polarons. 



5.3.4- Dynamical mean- field calculations Polaron formation has been extensively 
studied in the dynamical mean- field theory (DMFT) [73], which becomes exact for 
infinite degeneracy. Cappelluti and Ciuchi (149] developed a method for analytically 
solving the Holstein-t-J model in this limit. They found that antiferromagnetic 
correlations and the electron-phonon interaction mutually reinforce each other. 
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Figure 10. X/Xq as a function of J/t for ui p h/t — 0.1 and g/t = 0.1. The result 
of the SCBA (full line) is compared with the contribution from the coherent part of 
the Green's function (dashed) and the analytical formula ([30]) (after Gunnarsson and 
Rosch [T48] ). 

Cappelluti et al. |150j studied the relation between polaron formation and structures in 
the optical conductivity. 

The Holstein-Hubbard model in the paramagnetic phase has been studied 
extensively using the DMFT method (P-DMFT) pSH E521 [153]. It was shown that 
in the P-DMFT the electron-phonon interaction is strongly suppressed by the Coulomb 
interaction, in the sense that polaron formation is suppressed and states close to the 
Fermi energy are not strongly influenced by the EPI. The main effect was found to 
be a renormalization of U. As shown in Sec. 15.51 allowing for antiferromagnetic (AF) 
correlations by using an AF-DMFT method greatly increases the importance of the 
electron-phonon interaction. 

5.4- Differences between different phonons 

Sakai et al. |139j considered coupling via hopping integrals (in a three-band model) 
to breathing and buckling phonons involving in-plane and out-of-plane movements, 
respectively, of O atom in the Cu02 plane. By transforming to a t-J model, they 
obtained an on-site electron-phonon coupling 

H B b = uj ph b\ s b lS + g + b [s)( n i T n i+s ), (31) 

i,6 i,5 

where i labels the Cu sites and 5 = x,y differentiates the two bond directions. The minus 
sign between the two occupation numbers refers to a breathing phonon and the plus 
sign to a buckling phonon. For a breathing phonon, an O atom between two Cu atoms 
moves towards one Cu atom and thereby away from the other Cu atom, influencing 
Zhang-Rice singlets centered around the two Cu sites in opposite ways, as described by 
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Figure 11. Kinetic energy E^ n per hole as a function of coupling g for a CusOi6 
cluster with w p h = 0.2. All energies are in units of the hopping integral t. Open (solid) 
symbols correspond to the breathing (buckling) phonon. Solid (dashed) lines refer to 
two (one) doped hole. The figure shows how the breathing phonons reduce hopping 
more efficiently than buckling phonons (after Sakai et al. |139j ). 

the minus sign in Eq. ( I3~TI) . A buckling mode leads to a movement perpendicular to the 
CuC>2 plane of a O atom, which influences the two neighboring Cu sites in the same 
way, giving a plus sign in Eq. ( 13~TT) . The coupling to buckling phonons vanishes to linear 
order for a single perfect Cu0 2 plane, but it is finite for a plane with a static buckling. 

Figure [TT] shows the kinetic energy £<kin &s a function of the coupling constant g for 
breathing and buckling phonons. The figure illustrates that as g is increased, \E^ SX \ is 
more rapidly reduced for breathing than buckling phonons. If a (local) breathing phonon 
is excited, the potential is lowered on one site but increased on a neighboring site. This 
tends to strongly inhibit hopping and reduce l-E^inl- Exciting a buckling phonon, on the 
other hand, lowers the potential on both sites involved, and this inhibits hopping less. 
This illustrates the importance of what type of phonons the electrons couple to. 

Fehske et al. [M] observed similar effects when comparing the one-dimensional 
Holstein model (HM) to an extended Holstein model (EHM) with a Frohlich type long- 
ranged electron-phonon coupling decaying as l/d 3 for large d, where d is the distance 
between the electron and a phonon. Figure H21 shows the optical conductivity a Tes (u>) 
for the EH and EHM as well as the sum rule S reg (uj) = duJ a ieg (u'), where the Drude 
peak has been removed from <r reg (u;). The optical conductivity is determined by the 
current-current correlation function, where the current operator for the present models 
corresponds to the transfer of an electron to a neighboring site. In the strong coupling 
limit, an electron in the HM has a large binding energy, e p , essentially due to phonons on 
the same site. If the electron is moved to a neighboring site by the current operator, it 
looses the energy e p and it leaves behind excited phonons with the energy e p . Therefore, 
ct{uj) tends to have a broad peak centered at 2e p , which is at uj = 11 in Fig. [12] [64J. Due 
to the long-range of nature of the phonons in the EHM, phonons relatively far away 
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Figure 12. Optical conductivity er rcg (w) (thin lines) and sum rule S rcs (oj) (thick lines) 
for the Id Holstein (HM) and extended Holstein (EHM) models in the strong coupling 
limit (after Fehske et al. [64]). 



from the electron are excited in the EHM, even in the strong-coupling limit. Moving 
an electron to a neighboring site then costs much less energy, and cr(uj) has a peak at a 
much smaller energy [M]- Coupling to breathing phonons should instead shift the peak 
in a(u) to higher energies than in the Holstein model for a given A. 

There is a substantial coupling to apical oxygen phonons [HI [55], in particular to 
modes with a small q- vector parallel (qn) to the Cu02 plane. These phonons should have 
a similar effect as the phonons in the EHM model above, while modes with |q||| ~ n/a 
should be more similar to a breathing phonon. 

5. 5. Effects of Coulomb interaction 

Most of the work above used the Holstein- 1- J or large U Holstein-Hubbard model and 
focused on the effects of antiferromagnetic (AF) correlations. Here we focus on the 
Holstein-Hubbard model, which allows a continues increase of U and also to study other 
effects of U. 

In Sec. I5.3.4[ it was found that treating the paramagnetic state in DMFT 
(P-DMFT), leads to a strongly suppressed EPI. Due to the nature of P-DMFT, 
antiferromagnetic (AF) correlations are suppressed, which have been found to be 
important for the EPI (see Sec. 15. 3j) . This suggests that it would be interesting to 
apply an AF-DMFT method, where an AF state is allowed. A second reason for this 
is that the half-filled Holstein-Hubbard model must be an insulator for large U. In 
the P-DMFT this can only happen via Z —>■ 0. From Eq. ( 1301) it follows that this 
strongly suppresses the EPI, at least in the weak-coupling limit. In the AF-DMFT, on 
the other hand, it is possible to have an insulating state with Z > 0. A third reason 
is to notice that P-DMFT is equivalent to solving an Anderson impurity model (with a 
self-consistent host). We consider the electron Green's function, describing, for instance, 
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Figure 13. Z as a function of A for different U and for u>q — 0.01251^, where W is 
the band width. The figure shows how the Coulomb interaction moderately suppresses 
polaron formation [Z — > 0) (after Sangiovanni et al. 65J). 



the removal of an electron in photoemission. We focus on the corresponding final states 
close to the Fermi energy, i.e., in the Kondo resonance. These states have essentially 
the same occupancy of the local level as the ground-state, since the electron removed 
in the photoemission process is replaced by an electron hopping in from the host [154] . 
Actually, in the limit of infinite orbital degeneracy and an infinite U, the occupancy of 
the local level is unchanged [155] . Seen from the phonons, coupling to the net charge of 
the local level, a photoemission process corresponding to the Kondo peak then leads to 
no change. As a result, Holstein phonons have only an indirect influence on these states 
due to a renormalization of the parameters [1541 1151] . In the Holstein-Hubbard model, 
however, an electron filling the hole created in photoemission comes from another 3d 
level, which also couples to phonons, and in general this may be expected to influence 
the spectrum also close to the Fermi energy. 

The AF-DMFT method has been applied to the Holstein-Hubbard model on a Bethe 
lattice using exact diagonalization for solving the impurity problem [65]. The results 
for Z are shown in Fig. [13j The results for U — shows how Z is reduced from Z — 1 
for A = to Z « for A = A c = 0.33. We use this as the criterion for (small) polaron 
formation. For A = 0, an increase of U leads to a decrease of Z. However, Z deceases 
more slowly with A for a finite U, and polaron formation happens at a somewhat larger 
value A c . Thus the Coulomb interaction moderately suppresses polaron formation, at 
least in AF-DMFT. This is also shown in Fig. [TH 

It is important to notice that in the half-filled large U case, there is no polaron 
formation in the ground-state, since U suppresses charge fluctuations. The Green's 
function describes the final state after an electron has been removed (in, e.g., 
photoemission) and Z — > describes how the corresponding hole localizes due to 
polaronic effects. This is different from the Holstein model, where polaron formation 



Interplay between electron-phonon and Coulomb interactions in cuprates 



26 







U=0 U large 
Half Half (AF) 



1 



U large 
Half (F) 




U 



AF 



Figure 14. Critical value A c for polaron formation (Z — ► 0) in the half-filled 
Holstein model (U — 0) as well as in the large U Holstein-Hubbard model in the 
antifcrromagnetic (AF) or ferromagnetic (F) states. The figure illustrates that AF 
correlations help the EPI (smaller A needed for polaron formation) but that the net 
effect of U is a moderate suppression of the EPI. 



means the formation of polarons also in the ground-state. 

Macridin et al. [156] performed a dynamical cluster calculation (DCA) for the 
Holstein-Hubbard model using a 2 x 2 cluster. Using a different definition of polaron 
formation, they found A c > 0.5, similar to the result in Fig. [131 They emphasized the 
synergistic cooperation between the EPI and AF correlations. As a result they found 
that the AF transition temperature at finite doping is enhanced by the EPI. Macridin 
et al. [156J and Fu et al. [157j pointed out that the EPI can contribute to a charge 
density modulation seen experimentally. 

To see the effects of the antiferromagnetic correlations we compare with polaron 
formation in the ferromagnetic state (F in Fig. [Till , where antiferromagnetic correlations 
are completely suppressed. We then find that A c is very large. When the 
antiferromagnetic correlations are reintroduced, A c is strongly reduced (see Fig. HI]) , 
meaning that the electron-phonon coupling becomes more efficient, as expected. 

Although the antiferromagnetic effects, caused by U, strongly reduce A c , the net 
effect of U is still an increase of A c , due to other effects of U. We may then ask what 
this is due to. The Green's function of the ferromagnetic half-filled state, describes the 
creation of a hole in an otherwise filled spin up band. In the absence of phonons, this 
hole could move completely freely, and if the system has electron-hole symmetry, this 
state is equivalent to a Holstein model with a single electron at the bottom of the band. 
The half-full Holstein model and the single electron Holstein model were compared in 
Sec. 13.31 It was shown that A c is much larger in the single electron case, because the 
absolute value of the hopping energy per electron is much larger. Counterintuitively, A c 
is therefore larger for the ferromagnetic Holstein model than for the half-filled Holstein 
model because the hopping energy of the hole to be localized in the ferromagnetic case 
is larger, although the total hopping energy is strongly suppressed. 

The AF-DMFT calculation can also be applied to the doped system. Figure [Pol 
shows Z as a function of A for different fillings n. As the filling is reduced (doping 
increases) the critical A c for polaron formation increases. The reason is that AF 
correlations decrease with increasing doping, which reduces the effects of the EPI. This 
is consistent with the experimental observation [7] that polaron formation is gradually 



Interplay between electron-phonon and Coulomb interactions in cuprates 



27 




Figure 15. Z as a function of A for different magnetic moments m and associated 
fillings n for U = 1.75W and ujq = 0.01251^. The figure illustrates how the critical 
A c is increased as the filling is reduced (doping is increased) due to a reduction of 
antiferromagnetic correlations (after Sangiovanni et al. [65j). 




Figure 16. Configurations involved in the formation of the Zhang-Rice singlet. 



suppressed as the system is doped. In addition to the effect discussed here, the calculated 
strong coupling to apical oxygen phonons [55] becomes more efficiently screened as the 
system is doped, also reducing the tendency to polaron formation. 

5.6. Coupling constants 

Above we have discussed extensively how the EPI is influenced by the Coulomb 
interaction and antiferromagnetic correlations by comparing the Holstein-t-J or 
Holstein-Hubbard models with the Holstein model, assuming that the coupling constants 
remain the same as U is increased. As discussed below, however, the coupling constants 
themselves can change in an essential way. 

5.6.1. Breathing phonons Starting from a three-band model [35J, a t-J model with 
phonons can be derived [H]. A similar derivation can be made for the case when U — 
by projecting out the oxygen 2p-states. The EPI is different in the two cases, due to 
the formation of a Zhang-Rice singlet in the large U case. The coupling to a singlet 
involves extra prefactors. The main reason for the difference, however, is illustrated 
in Figure [161 The undoped ground state is nominally a d 1 state in the hole picture. 
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Doping adds a 2p-hole, nominally leading to a d 1 p 1 configuration. This configuration 
couples to two configurations, d°p 2 and d 2 p°. The corresponding projection (Lowdin 
downfolding) for the U = case only results in the coupling to one configuration [22J. 
The two coupling possibilities in the large U case lead to an enhancement in the coupling 
constant squared of the order of three. In addition the coupling is screened in the U = 
case, which reduces the coupling by an additional factor of about two [22] . This suggests 
that many-body effects are crucial for the coupling to the breathing phonons. 

5.6.2. Apical oxygen phonons Neutron scattering experiments [5] and calculations [55] 
suggest that there is a strong coupling to apical oxygen phonons. This strong coupling 
was found [55] to depend crucially on the poor screening in these systems, in particular 
for the undoped system, being an insulator. This was studied further by Meevasana 
et al. |104[ 1158] . They extracted the effective EPI for optimally doped Bi2201 from 
the structures of the PES spectrum in the nodal direction and compared this with the 
experimental loss function in the c-direction, obtained from optical measurements. A 
rather good agreement was found and it was concluded that an essential part of the 
coupling is due to c-axis O phonons. For the overdoped sample, they then argued that 
phonons below an energy of the order of 60 meV are screened and do not essentially 
contribute to the EPI. The result is then a reduction of the coupling and a stronger 
emphasis on high-lying phonons, in agreement with the EPI extracted from experiments 
for overdoped Bi2201. 

6. Phonon spectral function 

6.1. Phonon softening 

The softening of the (half-)breathing phonon has been studied by several groups 
[3S1 E21 HQ1 BH H21 S31 HI]- Von Szczepanski and Becker [38] derived a t-J model 
with phonons starting from a three-band model. They calculated the density response 
function for the electronic system using exact diagonalization for a small cluster and 
from this the phonon self-energy and softening. Khalliulin and Horsch [39j HDJ UH H2] 
calculated the density response of the t-J model using both slave bosons [39] and slave 
fermions [JT] , from which they deduced the phonon self-energy. Their calculated relative 
softening for the half-breathing [q = (it, 0)] and breathing [q = (tt, tt)] phonons are 
compared with experiment in Fig. [171 The theory correctly predicts that for intermediate 
dopings the half-breathing [q = (tt, 0)] phonon is softened more than the breathing 
phonon, although the coupling constant is larger for the breathing phonon [Eq. (jHJ)]. 
For large dopings (5 ~ 0.3) the softenings of the two modes become comparable. This 
was predicted before the experiment had been done. 

Rosch and Gunnarsson [H] derived a. t-J model with phonons starting from a three- 
band model. Using input from ab initio calculations they obtained the electron-phonon 
coupling. The t-J model was solved using exact diagonalization, including the phonons 
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Figure 17. Relative softening of the half-breathing (tt, 0) and breathing phonons 
("7r, 7r) as a function of doping 5. The left figure shows theoretical results of Horsch 
and Khaliullin [42] together with experimental results known at the time of the 
calculations and the right figure shows experimental results of Pintschovius et al. [6] 
(after Pintschovius [6]). 




Figure 18. Phonon dispersion in the (1,0) and (1,1) directions. Experimental results 
(dotted line) for x — and x = 0.15 are shown. Theoretical results (full curve) 
for x — 0.125 show the calculated softening from the experimental x = results. 
The average over boundary conditions is shown and the bars show the spread due 
to different boundary conditions. There is a strong softening in the (1,0) direction, 
while the softening in the (1,1) direction is weaker at this doping (after Rosch and 
Gunnarsson [4"4"]). 



in the calculation, and the phonon spectral function was calculated. Their calculated 
dispersion for S = 0.125 is compared with experimental results in Fig [181 The dotted 
curves show experimental results and the full curves the calculated softening for a 4 x 4 
cluster. The bars show the spread of the results due to different (periodic, antiperiodic 
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Figure 19. Phonon spectral function for different values of A. The bare phonon 
frequency is w p h = 0.1D and a Lorentzian broadening with the full width half maximum 
of 0.041? has been introduced, where D = W/2 is half the band width. The figures 
a-c show how the phonon softening at half-filling is dramatically suppressed by U and 
figure d that the softening increases with doping S (after Sangiovanni et al. [65j). 

or mixed) boundary conditions. The figure illustrates that the softening is large for the 
(1,0) than the (1,1) direction for this doping. The softening in the (1,0) direction is large 
for |q| > 7r/(2a). The softening essentially follows the coupling strength ~ sm 2 (q x a/2) 
but is larger at q x = n/{2a) than would be expected from this argument. 

Sangiovanni et al. [65] calculated the phonon softening within the AF-DMFT theory 
for the Holstein-Hubbard model. Figure fl~9a shows that there is a large softening for 
the undoped system if U is small. Figure fl~9b-c show how the softening is much smaller 
for larger values of U, although larger values of A were used in Fig. [T9~b-c. In particular, 
in Fig. [T9~b almost no softening is observed. Finally, Fig. IT9H shows results for a doped 
system. As the doping is increased, the softening increases, although U is large. This 
illustrates the sum rule in Eq. (j23|) for phonons. 

Falter et al. [89, 90J have developed a model for the charge response of cuprates, 
including the ionic nature of the system. Based on this model they have calculated 
the softening of several phonons. In particular, they predicted the softening of the Of 
phonon (see Fig. CP [891 EQ] before it was observed experimentally. For Lai.ssSro.isCuC^ 
the softening of the Of phonon is about 30% and its width is about 4 THz=17 meV 
[5]. This large softening and width suggests a strong coupling to doped holes. 

6.2. Phonon width 

Khaliullin and Horsch [40J have calculated the width of the (half-)breathing phonon in 
the t-J model and found that it is very broad for q = (71", 0) due to the coupling to a 
low-lying collective mode in the density response function, while the q = (n, tt) phonon 
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is narrower for optimum doping. This is in good agreement with experiment. 

It is interesting that LDA calculations [20] predict the frequency of the half- 
breathing phonon of YBa2Cu307 quite accurately, while the theoretical width is an order 
of magnitude smaller than the width measured for La 2 _ a: Sr a .Cu04 [5]. Since the shift 
(real part) and width (imaginary part) are both determined by the phonon self-energy, 
and since the real and imaginary parts are related via the Kramers-Kronig relation, this 
is a surprising result. This was addressed [22J by projecting the three-band model onto 
a one-band model using either the Hartree-Fock (HF) approximation or by including 
many-body effects using the Zhang- Rice method [36] to obtain the t-J model. The 
HF approximation shows similarities to the LDA approach. As discussed in Sec. 15.6.11 
many-body effects enhance the EPI coupling constants in the derivation of a t-J model. 
On the other hand, based on a sum rule in Sec. 15.1.11 many-body effects lead to a 
strongly doping dependent suppression of the EPI for the phonon self-energy. These 
two effects were shown to roughly cancel for 5 ~ 0.15 for the phonon softening [22]. Due 
to many-body effects, however, there are particularly many low-lying excitations which 
couple strongly to the half-breathing phonon and lead to a large width for this phonon 
[221 ESI Hd SD E2]. This effect is not present in the HF approximation [22]- 

7. Electron spectral function 

7.1. Polaronic behaviour in the undoped system 

7.1.1. Quasiparticle weight in the absence of phonons In Sec. 14.31 we presented the 
arguments of Shen et al. [7\ that undoped cuprates show polarons. An essential part 
of the argument was that the quasiparticle cannot be seen experimentally, because its 
weight Z is so strongly reduced by the interaction with the phonons. However, one 
could also imagine that coupling to spin fluctuations alone could have this effect, and 
there would then be no need to invoke phonons. We therefore first discuss whether or 
not Z is finite for the undoped t-J and Hubbard models. 

Some approximate calculations for the t-J and Hubbard models gave Z = for 
5 = |1591 1160] . Similarly, DMFT calculations [73J for the paramagnetic state of the 
Hubbard model gave a very small Z for a small 5 and a large U [1521 1153j . although 
later calculations including antiferromagnetic correlations using the AF-DMFT gave a 
substantial Z [161] . It is interesting that in the self-consistent Born approximation the 
carrier couples to magnons whose energies go to zero for |q| — > 0. There is, however, no 
infrared singularity, since the coupling also goes to zero, and the quasiparticle weight 
converges to a finite number as the system size goes to infinity [146J. 

Exact diagonalization calculations for undoped t-J clusters with 16, 18, 20, 26 
and 32 sites obtained finite values for Z, and there was no sign of Z going to zero with 
increasing cluster size [24 11162] . Brunner et al. [163] studied a 24 x 24 t-J cluster using a 
loop algorithm and extrapolated the results to the thermodynamic limit. Mishchenko et 
al. [164] calculated the Green's function Gk{T) for imaginary times r at h = (tt/2,tt/2) 
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Figure 20. ARPES spectra for the undoped system at T = for different k normalized 
to the height of the phonon side band and as a function of the binding energy (BE) 
(after Rosch and Gunnarsson [55]). 

for a 32 x 32 t-J cluster. By using a continuous-time worm algorithm [165] they could 
eliminate any systematic errors. From the r dependence of G*;(t) and from analytical 
continuation |166j they obtained Z. Both groups [1631 [L64J found values of Z similar 
to what had been obtained by exact diagonalization for small clusters. These results 
strongly suggest that Z stays finite for 5 = 0, and that therefore the lack of a visible 
quasiparticle in ARPES for undoped cuprates is due to the electron-phonon coupling. 

7.1.2. Effects of phonons In view of the results in the previous section, it is then natural 
to ask if the electron-phonon coupling is strong enough to give polaronic behaviour. 
This was studied for La2Cu04 [55] within the t-J model together with a shell model 
[75] for describing the phonon eigenvectors. From the eigenvectors one can calculate the 
electrostatic potential acting on a Zhang- Rice singlet due to the excitation of a phonon, 
which provides an essential part of the coupling. In addition, the phonons modulate the 
Cu-0 hopping integrals and the energy difference between the Cu and O levels in the 
three-band model used to derive the t-J model. This leads to an additional coupling 
mechanism. Defining the coupling as 

X = 2J2M 2 /(^N), (32) 

A = 1.2 was obtained |55j. This is well above the values A c = 0.4 for the Holstein-t-J 
model [63] and A c = 0.55 for the Holstein-Hubbard model [65] giving small polarons. 

We find that the dominating coupling for La2Cu04 is due to the (half-)breathing 
phonons, several apical oxygen phonons and some low-lying modes involving mainly La 
and Cu atoms. This is supported by inelastic neutron scattering experiments, showing 
a large broadening and softening under doping for the (half-)breathing and Of apical 
oxygen phonons. 
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Figure 21. Dispersion of the phonon side band (filled circles) and the quasiparticle 
(filled squares) for A = 0.46 as well as the dispersion of a broad peak (open circles) and 
the quasiparticle (open squares) for A = 0.4. The two curves show the dispersion of a 
pure t-J model ( |169|fT45] ) with the energy zeros defined appropriately for comparison 
with A = 0.46 (full line) and A = 0.40 (dotted line) (after Mishchenko and Nagaosa 
[63]). 



The resulting t-J model with phonons was solved using a method based on a 
statistical sampling of the phonons and exact diagonalization [551 H67[ 1154] . The results 
are shown in Fig. [20l The binding energy is measured in relation to the quasiparticle, 
which is too small to be seen. The peaks at about -1.1 to -1.3 eV are phonon side bands 
due to many unresolved phonon satellites. The width (0.5 eV) of the side band is close to 
the experimental result (0.47 eV) [55j, and also the T-dependence is in agreement with 
experimental observations [551 1 168] . The binding energy of the phonon side band is too 
large compared with experiment, which may indicate that the calculation overestimated 
the coupling strength. Reducing all coupling constants by a factor 0.8 (giving A = 0.75) 
gives the width 0.4 eV and the binding energy 0.6 eV, in reasonable agreement with the 
experimental results 0.47 eV and 0.5 eV, respectively. These results suggest that the 
EPI is sufficiently strong to give the polaronic behaviour seen experimentally. 

As discussed in Sec. 14.31 the dispersion of the main peak in ARPES spectra of 
undoped cuprates can be well described by the dispersion of the quasiparticle in an 
extended t-J model with up to third nearest neighbor hopping [26]. However, this 
model cannot explain the large width of the peaks. Shen et al. [7] therefore proposed 
that the peaks are actually not quasiparticles but phonon side bands. This then raises 
the question of why the phonon side bands should disperse like the quasiparticles in a 
model without phonons. Mishchenko and Nagaosa [63j performed diagrammatic Monte- 
Carlo calculations for a t-J model with and without phonons. Their results are shown 
in Fig. [2TJ For A = 0.46 there is a well developed phonon side band (filled circles) which 
is found to very closely follow the dispersion of the quasiparticle peak for A = (full 
line), strongly supporting the interpretation of Shen et al. [7J. A simple explanation 
of the results of Mishchenko and Nagaosa [63] has been given [167J. As expected, the 
weak quasiparticle (filled squares) for A = 0.46 shows almost no dispersion. For A = 0.4 
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the phonon side band is less well developed. 

7.2. Differences between phonons and spin fluctuations 

The t-J model with one hole can be replaced by a model where the hole couples to 
magnons, treated as bosons [143[ I144[ 11451 11461 1147j . In this approximation, phonons 
and spin fluctuations are treated on the same footing and can be directly compared. In 
particular, we can compare the couplings as defined by Eq. (1321) . The coupling to the 
magnons is then X M — t/(2J) = 1.67 for J/t = 0.3 while the coupling to phonons is 
only A = 1.2 or 0.75 according to the estimates in Sec. 17.1.21 It is then interesting to 
ask why not spin fluctuations alone can drive polaron formation. 

Mishchenko and Nagaosa [63] and as well as Ciuchi et al. [671 H49j pointed out 
that in a diagrammatic description, polaron formation requires the inclusion of diagrams 
where lines describing phonon Green's functions cross. Liu and Manousakis |I 147] showed 
that due to symmetry reasons, whole classes of such diagrams are identically zero for 
coupling to magnons. This should be an essential reason why magnons alone cannot 
lead to polarons in the sense that the quasiparticle weight goes exponentially to zero as 
is the case for coupling to phonons. Alternatively, one can notice that polaron formation 
for Holstein phonons involves the excitation of many phonons on the same site as the 
hole. The magnons are due to flipping spins with s = 1/2, and locally it is only possible 
to flip such a spin once. This is in contrast to local phonons, which are true bosons and 
can be excited infinitely many times. 

7.3. Kinks 

Much of the interest in the EPI in the context of cuprates was triggered as Lanzara 
et al. [S\ emphasized the presence of kinks in the dispersion of the photoemission 
spectrum (see Fig. [5]). There have been a large number of theoretical studies of this 
effect [HOI QUI Q221 El Q251 EH! E21 QZZ]. An extensive theoretical study was 
performed by Sandvik et al. [59], who solved the Eliashberg equations approximately, 
considering coupling to Holstein, breathing and buckling phonons as well as to the 
resonance peak seen in neutron scattering. The q dependence of the coupling for these 
modes is assumed to be rather different. The calculated spectra, however, did not show 
large qualitative differences. It was concluded that from the q dependence alone, it 
might be hard to determine which mode causes the main coupling |59j . 

Devereaux et al. [52] performed similar calculations, focusing on the coupling to 
breathing (at u Br = 70 meV) and B l9 buckling (at u Blg = 36 meV) phonons and 
comparing with experimental results for Bi 2 Sr2Cao.92Yo.o8Cu208+<5. They argued that 
the Big phonon in particular couples to the antinodal point while the breathing phonon 
couples mainly to the nodal point. In the superconducting state, they found a structure 
in the antinodal direction at about us lg + A = 71 meV, where A = 35 meV is the gap, 
while in the nodal direction they found a structure at about = 70 meV. The effects 
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of the c-axis phonons on the kink has been studied by Meevasana et al. |1U41 1158] , 
as discussed in Sec. 15.6.21 

Norman et al. |105] noticed that the peak-dip-hump feature seen at the antinodal 
point below T c could be explained by a coupling to the resonance peak, since this peak 
also only appears below T c , has the right energy and would couple particularly strongly 
to the antinodal point. It was concluded that there may be a strong coupling to this 
resonance. This led to a substantial amount of work [170[ I172[ 11731 I174[ 11751 1 177] 
describing a kink in terms of the coupling to spin fluctuations and to the resonance 
peak. 

The kink is particularly pronounced in the anti-nodal direction in multilayer systems 
below T c . The resonance peak is observed under these conditions, and it may therefore 
contribute to the kink. For symmetry reasons, the Bi 9 buckling phonon has a substantial 
coupling for multilayer systems [52J. The pile up of density of states around the 
superconducting gap should enhance the contribution of the B ls phonon to the kink for 
T < T c [52J. This makes it hard to determine the relative importance of the resonance 
peak and the B l9 phonon for the kink in the anti-nodal direction in multilayer systems. 

Eschrig and Norman [173] observed that since the resonance peak is seen for bilayer 
systems in the odd channel it couples bonding (with respect to the two layers) to 
antibonding states. Since there is a large density of antibonding states close to the 
Fermi energy, the resonance peak should in particular influence the bonding state. They 
concluded |T73j that the spectra of Bi2212 close to the antinodal point, interpreted in 
terms of antibonding and bonding states with a peak-dip- hump structures |176j . could 
be described if a sharp mode was introduced in the odd channel but not in the even 
channel. An interesting question is what happens if there is a sharp mode in both 
channels, as expected for phonons in the CuO"2 plane. Borisenko et al. |178j applied 
the same idea to the spectra close to the nodal point. They found that the width of a 
bonding state is larger than for an antibonding state of the same energy (but different k), 
and concluded that this favors coupling in the odd channel. If, however, the broadening 
at zero binding energy (perhaps due to impurity scattering and other elastic effects) is 
subtracted, the conclusion is actually reversed. 

7.4- Isotope effect 

There have been a number of theoretical studies of the isotope effect [1 791 1180[ 11811 1182[ 
11831 1184] . Andergassen et al. |179] studied the phase diagram related to a quantum 
critical point (QCP) due to incommensurate charge ordering. Going beyond a mean- 
field theory by including fluctuations, they found that the QCP can be shifted by the 
isotope effect, leading to a general shift of the phase diagram. This shows up as a 
strong isotope effect in various properties |179j . Seibold and Grilli |181] argued that 
the correlation length for charge order fluctuations might have an isotope effect and 
found that this can lead to isotope effects in the PES spectra. Paci et al. |180] studied 
the Holstein model and demonstrated a strong isotope effect in, e.g., the effective mass 



Interplay betwe* 



36 




X X 



Figure 22. T c versus A for for a Hubbard model with Holstein (HH), buckling (HBC) 
or breathing (HBR) phonons at (a) 5% and (b) 15% doping for k> p h = 0.3t and U = 8t 
(after Macridin et al. [186]). 

when the EPI was strong enough to put the system close to (small) polaron formation. 
In a study of the Holstein-Hubbard model, Paci et al. [184] found that the competion 
between the Coulomb and electron-phonon interactions strongly influences the isotope 
effect. Mishchenko and Nagaosa [183J studied the undoped Holstein- 1- J model using 
a diagrammatic Monte-Carlo method. They found that the isotope effect can be large 
under certain circumstances. 

8. Superconductivity 

Honerkamp et al. [1571 1185] have studied the effects of phonons on the superconductivity 
transition using a weak-coupling functional renormalization group method. They 
studied Holstein, breathing, A lg and B lg out-of-plane oxygen buckling modes. They 
found that the Holstein and Ai g phonons are unfavorable for <i-wave superconductivity. 
These phonons add an on-site interaction, which reduces the Coulomb repulsion. In 
the weak-coupling limit studied by Honerkamp et al. [1571 1185] . this was found to 
be unfavourable for d-wave superconductivity. Of the phonons studied, only the B l9 
phonon, which gives no on-site attraction for q = 0, was found to be favorable for 
d-wave superconductivity. 

Macridin et al. [1561 1186j performed DCA calculations for Hubbard model with 
coupling to Holstein, breathing and out-of-plane oxygen buckling phonons for a 2 x 2 
cluster. The calculations neglected the EPI on the hopping integrals and thereby the k 
dependence of the coupling. The results for the transition temperature T c are shown in 
Fig. [22j The figure shows how T c drops with A for all the phonon couplings considered. 
It was actually found that all phonons considered enhance the pairing, and it was 
speculated that this is due to an enhancement of the AF susceptibility [186] . The 
EPI, however, reduces the quasiparticle strength Z, and it was concluded that this is 
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more important than the increase of the pairing, leading to a reduction of T c [186J. It 
would be interesting to also study the effects of coupling to hopping integrals for the 
buckling modes and the effects for a larger U. 

9. Summary 

We have reviewed evidence that the electron-phonon interaction (EPI) can substantially 
influence various properties of the high-T c cuprates. Some of the evidence is indirect and 
it can then be hard to determine which phonon modes are involved or to distinguish 
between phonons and other bosonlike excitations, such as the resonance peak. It is 
therefore of particular interest that inelastic neutron scattering shows that certain 
phonon modes are appreciably broadened and shifted when the system is doped. 
This is strong evidence that these modes have a substantial EPI. In particular, such 
effects have been found for apical oxygen phonons and the (half-)breathing oxygen 
bond-stretching phonons. Theoretical many-body calculations indeed find substantial 
couplings for these phonons. For multilayer systems similar effects are also seen for 
other phonons, in particular, the B l9 phonon. We have discussed theoretical treatments 
of phonon softening and broadening and illustrated that the experimental behavior can 
be described theoretically. 

Much of the interest in the EPI was triggered by the observation of kinks in the 
dispersion seen in photoemission, which can be interpreted in terms of coupling to 
phonons. It remains controversial, however, how much phonons and the resonance peak 
or other spin excitations contribute to these kinks. Theoretical treatments of the kink 
have been shown. 

There is strong evidence in favor of (small) polaron formation in undoped cuprates. 
While this polaron formation is helped by the antiferromagnetic correlations, we have 
argued that antiferromagnetic fluctuations alone could not lead to small polarons. This 
suggests a substantial EPI for the undoped cuprates. The EPI for LaCu04 has been 
calculated, obtaining substantial coupling to apical and breathing phonons, and it was 
shown that the experimental line shape of the photoemission spectra can be understood 
rather well. As the system is doped, the polarons disappear. Photoemission spectra, 
however, still show substantial weight in the energy range where the polaron related 
phonon sideband was observed, suggesting that the EPI is substantial also in the doped 
systems. 

An important feature of the cuprates is the great importance of the Coulomb 
repulsion. We have therefore treated the interplay between Coulomb and electron- 
phonon interactions extensively. Using sum rules, we showed that for weakly doped 
cuprates, the Coulomb repulsion strongly suppresses the phonon self-energy, while there 
is no corresponding strong suppression of the electron self-energy. For polaron formation 
as seen in photoemission (described by the electron self-energy), antiferromagnetic 
correlations, resulting from the Coulomb repulsion, greatly help the EPI. Nevertheless, 
due to other effects of the Coulomb interaction, there is a moderate suppression of 
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polaron formation, at least in the antiferromagnetic DMFT. 

The EPI is usually discussed for some model with EPI, assuming that the electron- 
phonon coupling constants are fixed. We find, however, that when deriving, for instance, 
Holstein-t-J or Holstein-Hubbard models, the coupling constants may be substantially 
enhanced by many-body effects. This is, for instance, the case for the (half-)breathing 
phonons. 
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